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The low momentum behavior ol the Landau gauge Gribov-Zwanziger action is investigated using 
the respective Dyson-Schwinger equations. Because of the mixing of the gluon and the auxiliary fields 
four scenarios can be distinguished for the infrared behavior. Two of them lead to inconsistencies and 
can be discarded. Another one corresponds to the case where the auxiliary fields behave exactly like 
the Faddeev-Popov ghosts and the same scaling relation as in standard Landau gauge, ka + 2k c = 0, 
is valid. Even the parameter k is found to be the same, 0.595. The mixed propagators, which appear, 
are suppressed in all loops, and their anomalous infrared exponent can also be determined. A fourth 
case provides an even stricter scaling relation that includes also the mixed propagators, but possesses 
the same qualitative feature, i. e. the propagators of the Faddeev-Popov ghost and the auxiliary 
fields are infrared enhanced and the mixed and the gluon propagators are infrared suppressed. In 
this case the system of equations to obtain the parameter k is non-linear in all variables. 

PACS numbers: 11.10.-z,03.70.+k,11.15.Tk 



I. INTRODUCTION 



In recent years different methods have been used to determine Green functions, especially propagators, in Yang-Mills 
theory. A prominent choice of gauge is the Landau gauge due to its manifest Lorentz covariance and the minimum 
number of interaction terms, which renders it very accessible for functional methods. Whereas the ultraviolet parts 
can be obtained by perturbation theory, a more refined, i. e. non-perturbative, treatment is necessary for the low and 
intermediate momentum regime. The asymptotic behavior at large distances is thereby of special interest, since the 
confinement scenarios of Kugo and Ojima [TJ [5] and Gribov and Zwanziger [2HB] can be tested: The ghost propagator 
should be infrared (IR) enhanced, whereas the gluon propagator vanishes at zero momentum. 

The original Gribov-Zwanziger scenario is based on an improved gauge fixing compared to the standard Faddeev- 
Popov theory [TJ, which restricts the integration in field configuration space to the hyperplane d^A^ = 0. However, 
there are still configurations left that are related by gauge transformations, so that the Faddeev-Popov gauge fixing is 
insufficient. Gribov proposed by his no-pole condition a way to restrict integration to what is nowadays known as the 
Gribov region [3 . Zwanziger generalized the no-pole condition to all orders leading to the horizon function [1] , which 
is a non-local object. However, it can be localized by the introduction of more fields [I]. An immediate consequence 
is that the perturbative gluon propagator changes its form due to a mixing at the level of two-point functions with 
the new fields: Instead of having a massless pole it vanishes at zero momentum. The quantization along Gribov- 
Zwanziger requires a new parameter, the Gribov parameter 7, whose value is fixed by a gap equation. Taking into 
account this equation, one can show that the Faddev-Popov ghost propagator gets IR enhanced like l/k 4 [31 |5J [5] . 
Thus the improved restriction of integration in field configuration space alters the propagators significantly. This 
is a semi-perturbative result, since the perturbative propagators are considered but including the non-perturbative 
information the Gribov horizon provides. One should note that the UV part of the theory is not affected by the 
horizon function: In the asymptotic limit of very high momenta the propagators retain their usual form and also the 
renormalization of the theory requires no additional independent renormalization constants [51 1^1 ITU], 

Of course it is interesting to employ a completely non-perturbative method to the investigation of propagators. Such 
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a framework is provided by Dyson-Schwinger (DSEs )|TT] and renormalization group equations (RGEs) |12l [T5] . The 
former are the equations of motion of Green functions and consist of an infinitely large tower of coupled equations. 
Therefore results in principle depend on the truncation. However, within the so-called scaling solution it is possible to 
analyze consistently the qualitative behavior of all Green functions in the IR. Again the gluon propagator in Landau 
gauge is found to be IR vanishing, going like (p 2 ) 2K_1 , whereas the ghost gets IR enhanced as (p 2 ) _K_1 [HI 1 1 5] . 
The parameter k is positive and can be evaluated as 0.595 |16[ I17|. The behavior found by Gribov and Zwanziger 
corresponds to k — 1. One caveat of the use of Dyson-Schwinger equations up to now may be that only the standard 
Faddeev-Popov action has been used and Gribov copies were not taken into account explicitly. However, there exists 
an argument by Zwanziger, why the Dyson-Schwinger equations should not change, when derived with a proper 
restriction of the path integral |17j . The key to this is that the determinant of the Faddeev-Popov operator, which 
appears in the path integral after the standard gauge fixing, vanishes at the first Gribov horizon. Thus, if one formally 
cuts the integration there, no additional surface terms are introduced and the DSEs are not altered. 

The scaling type solution described above was the first solution found for both Landau gauge propagators using 
DSEs [14] [15]. A few years ago another solution was obtained, called the decoupling solution, where the gluon 
propagator freezes to a constant value in the IR and the Faddeev-Popov ghost remains unenhanced |18ff22] . This 
solution is favored by most recent lattice calculations [23 - 28J and the refined Gribov-Zwanziger framework [181 129"] 
and can also be obtained by functional equations [191 I2T1 122] . The scaling solution on the other hand is found in lattice 
calculations at /? = [30| and in two dimensions |31| . by DSEs and RGEs using appropriate boundary conditions 
[HI [TBI I2TI 1221 l3"2H3~i] and by stochastic quantization methods [T7l [5SJ [35] • So both solutions are found within 
functional methods, but it seems hard to get the scaling solution on the lattice. In this respect there has recently 
been a proposal |37_ on how to implement different boundary conditions as used in functional equations also on the 
lattice, so that both the family of decoupling solutions and the scaling solution can be obtained. 

The use of a non-perturbative method might clarify the connection between the Gribov-Zwanziger and the standard 
Faddeev-Popov Lagrangians. There have also been speculations that the negligence of Gribov copies in DSEs might 
be the reason for the observed discrepancy with lattice calculations. In this article we take a first step in this direction 
by analyzing the IR properties of the Gribov-Zwanziger action with Dyson-Schwinger equations. The basic techniques 
for such an endeavor are already known [151 138H41) , but due to the mixing of the gluon field and the auxiliary fields 
some modifications are necessary. This mixing is indeed the main source of intricacies arising in this context: When 
the propagator of a field and its two-point function are no longer the inverse of each other, but there is a matrix 
relation involving several propagators/two-point functions, it is necessary to distinguish between their respective 
infrared exponents. We deal with the non-linear relation between the different dressing functions by splitting our 
analysis into four different cases, depending on which part is IR leading. We will show that two of these cases lead to 
inconsistencies so that only two possible solutions remain. In both of them the ghosts and the auxiliary fields have 
the same qualitative IR behavior. Their infrared exponent is related to that of the gluon propagator by the scaling 
relation ka + 2Kghosts = 0. The differences are in the details. The first solution has an independent mixed propagator. 
Since its contributions in DSEs are suppressed, the effective system of equations reduces to the same system as in 
standard Landau gauge. Consequently also the parameter k := n g hosts > is the same. The infrared exponent of 
the mixed two-point function is also calculated and shows that the mixed propagator is IR suppressed. The second 
solution provides also for the infrared exponent of the mixed propagator a scaling relation according to which it scales 
with the exponent k/2 — 1 and is thus IR suppressed. 

Zwanziger found in a recent article an independent IR enhancement of the bosonic auxiliary field, where its propa- 
gator goes like (p 2 )~ 1 ~ d / 2 02]. His ansatz differs from ours as he starts with the propagators, whereas we start with 
the two-point functions. Furthermore his solution is based on a cancelation in the determinant of the propagators 
that appears in the relation between the two-point functions and the propagators. 

After setting the scene by presenting the usual Gribov-Zwanziger Lagrangian rewritten into a convenient form in 
sec. [IT] we extend the known methods for power counting to the case of mixed propagators in sec. |III[ The analysis 
is then carried out in sec. |IV| and we present our conclusions in Sec. [V] The appendix contains some details on the 
calculations of sec. IIIII 



II. THE DYSON-SCHWINGER EQUATIONS FOR THE GRIBOV-ZWANZIGER ACTION 

The Lagrangian of Yang-Mills theory in Landau gauge is 

C FP = \f; v F; v + ^(d,A«f - c a M ab c\ (1) 

with the field strength tensor 

FZ» = %K ~ 9uA% - g r hc A\Al (2) 
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and the Faddeev-Popov operator 



M ab := -d„Df = -5 ab n - f^A^d^ (3) 



The fields A, c and c denote the gluon, ghost and anti-ghost fields. £ is the gauge fixing parameter that is set to 
zero for Landau gauge, so that the integration in the path integral is restricted to the hyperplane d^A^ — in field 
configuration space. This restriction is sufficient, when only small fluctuations around A = occur (as in perturbation 
theory), since then no gauge equivalent configurations are taken into account. However, going to smaller momenta 
away from the ultraviolet, at some point so-called Gribov copies appear that falsify any calculations. These copies are 
gauge copies but within the hyperplane fixed by our choice of gauge. By restricting the field configuration space to 
the Gribov region, defined by those configurations that have only positive eigenvalues for the Faddeev-Popov operator 
[3], the problem gets alleviated, but still copies remain j?3]. To single out a unique representative of a gauge orbit 
one can for example take those configurations along the gauge orbit that are in the global minimum with respect to 
their norm. On a lattice this is, however, a polynomially hard computational problem and therefore not feasible on 
large lattices. How a restriction of the integration region to unique field configurations could be implemented at the 
level of the Lagrangian is currently unknown. So we ignore these copies for now and will only talk about the Gribov 
region in the following, for which such a restriction can be done. 

To accomplish this restriction one can add a non-local term to the Lagrangian called horizon function [4 . This 
introduces a new variable 7, the Gribov parameter. It is not free, but fixed by a gap equation. However, the non-local 
term can be localized using two pairs of auxiliary fields [3]: The bosonic fields ip and (p and the fermionic ones u> and 
Q. Here we give the localized expression, since we will only work with this form of the Gribov-Zwanziger Lagrangian: 

Cgz = (pl c M ab ^ - Qf°M°» u Z + ! 2 g f abc A^ - ft) - j 4 d(N 2 - 1), (4) 

with the Gribov parameter 7 determined by 

(9 f abc A a M c - (ft)) = 2 7 2 d(N 2 - 1), (5) 

where N is the number of colors and d the dimension of space-time. 

For our calculations we will split the bosonic ghost fields into real and imaginary parts U and V, respectively, as 
done in ref. |42| : 

£ = £ FP + £' GZ , (6) 
■Cgz = £u + £v + £uv - uMu, (7) 

Cu = \u™M ab Ul c , (8) 
C v = \v? M ab V bc + 1 g 7 2 V2f abc A a X c , (9) 



C uv = \igf abc ufV bd d v Al, (10) 



where U and V are defined by 



<P=A=(U + iV), (p=±={U-iV). (11) 



When the Landau gauge condition d^A^ = is enforced, Cjjv vanishes and the only mixing on the level of two-point 
functions is between the gluon field A and the imaginary part of the bosonic ghost field V, whereas the [/-ghost does 
not mix. This splitting simplifies calculations, because we only have to deal with a two-by-two matrix instead of a 
three-by-three matrix for the mixing. Perturbative expressions for the propagators derived from this Lagrangian can 
be found in ref. |44| . 

Now we have three different types of fields (Faddeev-Popov ghosts c and c, fermionic auxiliary fields to and Q, real 
part of the bosonic auxiliary field U) that only interact with the gluon field via the Faddeev-Popov operator. To 
alleviate calculations for our purposes it is useful that we can combine all these fields into one type only. Since they 
appear quadratic in the action, they can be integrated out in the path integral: 



D[cc}e 5Mc = detM, (12) 
D[Qu]e aMu = (detM)'* N *- 1 \ (13) 
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field 


# degrees of freedom 


statistics 


c, c 


1 


fermionic 


oj, Q 


d(N 2 - 1) 


fermionic 


4> 


d(N 2 - 1) 


bosonic 


u 


d/2{N 2 - 1) 


bosonic 


V 


d/2(N 2 - 1) 


bosonic 


V, V 


d/2(N 2 - 1) + 1 


fermionic 



TABLE I: The numbers of degrees of freedom and the statistics of the Faddeev-Popov ghosts c and c, the original auxiliary 
fields uj, oj, <j> and cj>, the bosonic auxiliary fields U and V and the fermionic fields r\ and f\. 

J D[U}e~i UMU = {detM)-^-^. (14) 

The different exponents of the determinant of the Faddeev-Popov operator are due to the different numbers of degrees 
of freedom. For 7 = also the U-field can be integrated out and all determinants from auxiliary fields cancel, so that 
the original Faddeev-Popov Lagrangian is recovered. For the purpose of this article we can treat all these non-mixing 
fields as new fermionic fields rj and fj with the appropriate number of degrees of freedom, namely |(iV 2 — 1) + 1. The 
field V cannot be included due to its mixing with the gluon field and therefore the two fields V and 77 can have a 
different infrared behavior. An overview of the different fields is given in tab. [I] The final Lagrangian reads 




where the subscript index of the new ghost fields rj and fj runs from 1 to |(A^ 2 — 1) + 1. Note that for odd dimensions 
and even N this number is half-integer and therefore this transformation is not directly possible. However, we can 
consider only integer values and perform an analytic continuation to half-integer values, if necessary. Alternatively 
one keeps the three ghosts separated and will get the appropriate numerical factors in front of the diagrams. A further 
diagonalization of the Lagrangian would require a further splitting of the fields in color space due to the different 
number of color indices of the A- and U-fields. The new fields would mean a significant complication of the Lagrangian 
at the level of vertices. Therefore we continue with the expression given above. 

The derivation of the DSEs for this Lagrangian is done in the standard way [TT1 143] . i. e. by variation of the path 
integral w. r. t. a field and further field differentiations. An intricacy is the mixing of the A- and U-fields, which 
allows already at the perturbative level additional vertices, e. g. an AAU-vertex. For the 77-field, on the other hand, 
no such vertices are allowed. The number of diagrams in the DSEs increases drastically by the mixed propagator and 
the additional vertices. The task of their derivation was performed with the Mathematica package DoDSE |45| . The 
full two-point equations are given in figs, [I] [2] [3] and [4] 

III. POWER COUNTING 

The IR part of propagators can be determined without solving the complete DSEs numerically. This technique 
relies on the assumption that at low momenta all dressing functions obey a power law with a specific exponent, called 
the infrared exponent (IRE) [15]. For small external momentum the integrals in the propagator DSEs are dominated 
by small loop momenta, so that all expressions can be replaced by the corresponding IR expressions, i. e. power laws. 
The integrals can be evaluated analytically or - if we are only interested in the qualitative behavior - it suffices to count 
the powers of all momenta in the integral, since every loop momentum will transform into an external momentum 
upon integration due to the lack of other mass scales. As a result the IREs of different propagators are related by 
so-called scaling relations, thereby reducing the number of independent IREs. 

This approach was very successful in Landau gauge, where all IREs can be expressed with one parameter k |14|,I15|. 
Using a bare ghost-gluon vertex its value is k = (93 - vT201)/98 « 0.595 [H E]. Varying dressing functions 
that respect the fact that the vertex does not have an IR enhanced dressing change this value only slightly |16j . 
Furthermore, the dependence on the kinematic configuration was also found to be mild [46 . The first calculations 
took into account only the one-loop diagrams in the propagator DSEs. Later on it was shown that this IR truncation 
indeed is consistent with the whole tower of DSEs [38 and also RGEs j33]. Therefore one combines the two systems 
of functional equations, where the decisive difference is that the DSEs feature a bare vertex in every diagram, whereas 
the RGEs only contain dressed quantities [34 , 39J . The generalization to an arbitrary Lagrangian was provided in ref . 
|41j . We shortly summarize this method, since a great part of the forthcoming analysis is based on it. 
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FIG. 1: The DSE of the gluon two-point function. A propagator with the exponent —1 denotes the two-point function. This 
convention, strictly speaking being mathematically incorrect, is chosen for the purpose of diagrammatic representation only. 
The propagators are labeled by the respective fields. The indices i and j denote the first and second fields of the depicted 
two-point function DSE. 

Before doing so we have to introduce some notation. The IRE of propagators and two-point functions are given by 
8ij and Kij, respectively. This distinction will be of importance later, because the relation between the two can be 
non-trivial if the Lagrangian contains mixed two-point functions, while for the case of non-mixed two-point functions 
it is Ku — —Su- Higher vertex functions get the IRE k^-.... To distinguish different propagators and vertices we add 
the corresponding fields as indices, e. g. ka C c for the ghost-gluon vertex. Also we refer to vertices and propagators by 
monomials of their respective fields, e. g. the ghost-gluon vertex is denoted by ^Ice-vertex. If we talk about DSEs the 
order of the fields corresponds to the order in which the differentiations have been performed. For instance AT^-DSE 
means the DSE of the vertex with two V-legs and one A-\eg, where the first derivative was performed w. r. t. the 
gluon field. This is important in some cases, when the order of the differentiations is crucial. Finally, we note that 
we use Ki := Ku for non-mixed two-point functions, e. g. ka ■— ha A- 

In the case of non-mixing two-point functions we get the relations presented below. The calculations are done for a 
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FIG. 3: The DSEs of the A-V-mixed two-point function. The first is the AV-T>SE and the second one the VA-BSE. 
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FIG. 4: The DSE of the ^-field two-point function. 
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set of fields {<f>i\. Considering a general diagram v with external legs of the field <j>i, the formula for its IRE k v is 
1 



2 ^ m A+ Yl <...vUE^^) + E ('••< ...< • ',E ; < ' ! 0- 

i vertices, r>3 \ i / vertices, r>3 \ i / 



where <5j and K^...^ are the propagator and vertex IREs, is the number of bare/dressed ^ . . . (/^-vertices and 

kY"' lr gives the number of 0,-legs of a (f)^ . . . 0j r -vertex. Furthermore it was found that DSEs and RGEs yield bounds 
on the last two terms, namely 



'"*"**> 0, (17) 



which exists for every possible vertex with the IRE Kj 1 ...i r , and 

~£*r-H>o, as) 

which only relates to the vertices appearing in the Lagrangian. Details on these two equations are given in Appendix 



A 2 Thus the last two terms in eq. (16 1 can be dropped to get the lowest possible IRE of the diagram v: 



3v,max — 2 ^ ^ Tflj&i ■ (1^) 

i 

It depends only on the number of external legs and is therefore the same for all diagrams of a DSE. Using this 
expression in the analysis of leading diagrams in the propagator DSEs, one arrives at the following inequality: 

\< E'^-; ' •"• ( 2 °) 

i 

where the hat indicates that this equation is only valid for the leading diagram(s). In every diagram of a DSE is 
exactly one bare vertex, so that there is only one n b = 1, while the others are zero. Since this inequality is just the 



opposite of eq. (18 1, it has to be saturated, i. e. 

1 



An immediate consequence is 



,»;'..,. E ! V--; ' (2i) 



= 0, (22) 



i. e. the vertex corresponding to n| i does not scale in the IR. This is a direct consequence of the Fischer-Pawlowski 
consistency condition namely that there is one bare vertex in the DSEs and none in the RGEs, but both equations 
should yield the same result. In this analysis it was not specified which diagram was assumed to be leading. This 



can be done now, i. e. one tries all possible realizations of eq. (21 1. Since there are only as many possibilities as 
bare vertices in the Lagrangian this is a manageable task. In other words, one can read off possible scaling relations 
directly from the Lagrangian and take only those that yield non-trivial results [41J. Let us illustrate this approach 
in standard Landau gauge: We only have three bare vertices: the ghost-gluon vertex, the three-gluon vertex and 
the four-gluon vertex. The two gluonic vertices yield the same inequality, which even simplifies our task. The two 
possibilities that remain arc: 

5 a = (from gluonic vertices) , (23) 
or 

8a + 25 c = (from ghost-gluon vertex) . (24) 

The first possibility corresponds to the perturbative (trivial) solution, i. e. all Green functions are bare, whereas the 
second solution yields the known scaling relation. As this simple example shows, the method presented in ref. |4"T] 
provides an easy method to derive a scaling relation. Its use is especially convenient for more complicated Lagrangians 
like Yang-Mills theory in the maximally Abelian gauge |41j . 
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The task is now to generalize these results to a Lagrangian with mixed propagators. An important consequence 
of this mixing is that a propagator is not the inverse of the two-point function, but there exists a matrix relation 
between propagators and two-point functions T^, 

D 4>4> = (p^)- 1 , (25) 

where cj> can be any field. This is the reason why we assign different IREs to propagators and two-point functions, 
namely 5ij and Kg, respectively, and lies at the heart of the complication of the method. 



Whereas eqs. ( 17 1 and ( 18 1 remain valid, as exemplified in Appendix |A| eq . ( 16 1 undergoes some modifications, since 
the topological relations used in its derivation have changed; see eqs. ( |A2| and (A3 1. The details of the calculation 
are given in App. [A] From now on we work with general dimension d. As a result we obtain for the lower bound on 
the IRE of a diagram 

2 ^2 m A + ^n AV (2S A v - 5 a - 5 V ). (26) 

i 

The new term is the last one, by which the formula depends on the number of internal mixed propagators. Also 




eq. (20 I gets modified: 



j V 3 / \ dressed 

vertices 



(27) 



where k A 1 y lr indicates the number of times a mixed propagator is contained in the diagram that determines the IRE 
of the vertex (f>i t . . . <f>i r . 

An issue we have to address is the role of the bare two-point functions appearing in the two-point DSEs. In the 
following analysis we ignore them. For the gluon propagator this makes no difference at all in the analysis, while for 
the other propagators this is of vital importance. The reason is that including the bare two-point function would 
restrict the IRE of the two-point function to negative values. However, for the Faddeev-Popov ghost propagator, the 
fermionic auxiliary field and the real part of the bosonic auxiliary field it is well justified to neglect it in an IR analysis 
due to the horizon condition, which can be used to show that the ghost two-point function has to vanish at zero 
momentum and thus the bare two-point function has to be canceled [5] 1421 147| . This is the way the horizon condition 
enters in our calculations. A similar mechanism for the mixed propagator is not known. We comment on this and the 
cancelation of the imaginary part of the bosonic auxiliary field V at the ends of subsec. |IVB| and |IVC| 

IV. PROPAGATORS AND VERTEX FUNCTIONS IN THE INFRARED 

Normally one does not have to distinguish between different tensors for power counting, since automatically the 
most divergent tensor is taken into account. However, when the relation between the propagators and two-point 
functions is not trivial as here, more care has to be taken. To keep the number of tensors manageable we use the 
following truncation of the tensor basis: We only take into account those tensors of two-point functions that appear 
in the Lagrangian. We will see that already in this minimal truncation a non-trivial V^-propagator arises, in the 
sense that there are two independent parts with different dressing functions. Thus they can behave differently and 
both of them have to be taken into account. 

Explicitly the matrix of the mixing two-point functions is given by 



(28) 



T AA,ac = 5 ac p 2 c ± (p 2 )p ^ + S ^- C \(p 2 )p^ (29) 



T vv,abcd = 5 ac S bd p 2 c v (p 2 ) 9flv , (30) 
T Av,cab = r Hp 2 c AV {p 2 )g^ (31) 
and that of the propagators by (£ = 0) 

/ r\AA r)Av\ 

D **=[ D VA D D vv)> (32) 
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T\AA.ab xa& 

" v " P 2 ^ci(p2)c v (p2)+2Nc\ v (p2y [66> 

j^VV.abcd _ 1 racsbd _ rabe rede }_ p 2c AV (p 2 ) 

p 2 Cy(p 2 )° 5 ^ ' ' ^^ c ± (p 2 )c ^ (p 2 )+2A r C ^(p2 )cv(p2) - ^ 



r-)AV,abc ■ rabe 

%I p 2 » v c\{p 2 )c v {p 2 ) + 2Nc 2 AV {p 2 V i "" 
For the fermionic ghost the standard relation is valid: 



Df d ab = {Tf^y 1 = -S^dJ^. (36) 



The functions are the dressing functions of the two-point functions with the IREs Kij defined as 

Cii(p 2 ) (37) 

Similar to the IREs we define da — di and cu — Ci . Because of the matrix inversion the determinant of the two-point 
function matrix appears, detC = c A (p 2 )cv(p 2 ) + 2N c AV (p 2 ), which forces us to distinguish between several cases: 

I: c 2 AV > cacv <-> ka + kv > 2kav 

II: CaCv > c 2 AV <H> 2K^y > KA + Ky 

III: c^y ~ c^cy f> + Ky = 2k j 4\/. no cancelations 

IV: c AV ~ c^cy •H- + = 2kavi cancelations 

If the IREs of the two expressions in the determinant are the same, it is possible that the IR leading parts cancel. 
Thus we have to take into account this case explicitly as case IV. 

We will now determine the form of the propagators for each of the four cases. 



A. Case I: ka + Ky > 2kav 



By assumption the dressing function cav dominates in the numerator. Therefore the propagators take the following 
form: 



D 



AA,ab %ab 



= s ab \p, 



c v {p 2 ) 



£)VV,abcd 

^ V 2 cv{p 2 ) 



p 2 ^2Nc 2 AV {p 2 Y 
1 1 



rac %bd _ rabe rcae p 



abe rede ; 



D 



AV.abc 



— _i f abc — r> 

~ 1 J „9 r i 



1 



N ' 



P 2 ^V2Ncav(p 2 ) 



(38) 

(39) 
(40) 



Note that both tensors contribute to the V ^-propagator, since they have the same IRE. The relations between the 
IREs of propagators and two-point functions are 



5a = k v - 2k av, 
<5y = — «y, 
Sav = —k-av- 



With these expressions we can determine the value of the additional term that appears in eqs. (26 1 and (27): 

2^Ay — &A ~ $V = — 2K^y — Ky + 2K^y + Ky — 0. 



(41) 
(42) 
(43) 



(44) 



Thus the two equations stay the same as in the case with no mixing propagators. The lower bound is then as usual, 



= u- 2 ) \ l -\H m A-\H m ^ 



(45) 
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and eq. (27) takes the form 

* + ij>m 3 - \(l-^{r-2) + ^ky-^5A >0. (46) 

It is now a proof of only a few lines to show that case I corresponds to an inconsistent assumption. In the V^-DSE. 
see fig. [3j we only have bare AW- vertices. Thus eq. (46) becomes 

KAV + \{&A + 6v)-\{6A + Mv)-\(i-2\>Q => 2k av >-k v + ^-2. (47) 
From the VF-loop of the AA-BSE, depicted in fig. [5] we get 

k a < -2kv + 7,-2 =>■ k a + n v < -kv + ^ - 2, (48) 

where kavv < from the AUU-DSE was used. Using these expressions in the defining assumption of case I, we find 
an unsatisfiable inequality: 

d „ „ d n 

—k v + - — 2 > k a + k v > 2k av > —ny + — — 2, 

> 0. (49) 

Thus case I can be disregarded except one potential caveat. For case I we cannot exclude that the mixed two-point 
function might not be determined by the loop integrals. Therefore we also discuss the case in which the tree-level 
term is the IR dominant part, which corresponds to n A y — —1. since the tree-level expression is constant. In the 



UU-DSE there is only one bare vertex, the AW- vertex. Thus eq. (46 1 yields 



2 

The yU-DSE itself leads from the AA-W-Xoa^, as given in fig. [5] to 



5 A + 2S V + % - 2 < 0. (50) 



2 

Together these two inequalities give 



S A + 2S V + - 2 > 0. (51) 



S A + 2S V + | - 2 = 0. (52) 

Inserting the expressions for the propagator IRE, we arrive at 

d d 
k v = - - 2 - 2n AV = -■ (53) 

The last inequality we need we get from the AA-DSE, namely from the UF-loop as depicted in fig. [5] 

k a < 25 v + ^ - 2 = -2k v + ~ - 2. (54) 

Now we insert these equations into the defining relation of case I: 

d „ 

-Ky + - — 2 > ka + Ky >2kav 

-2 > - 2. (55) 

We see that also the case of a constant two-point function leads to an inconsistency. 
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V V^, r A A-x 

a m > — a — . ■ — v—m > — v 

FIG. 5: Diagrams of the AA- and W-DSEs used explicitly in the text. 



B. Case II: ka + Kv < 2kav 



Here the dressing function combination c A cy dominates in the numerator: 

1 1 



D AA, ab = 5ab ~p __ (g6) 



n VV,abcd _ 1 ^ / A ac A bd n rabe rede p ^ c Av(P ) A . ^ j; , - - , 

p 2 C V {p 2 ) \ P C^(p 2 )Cy(p 2 )y p 2 Cy(p 2 > 



1 D \/2c A v(p 2 ) 



Only the first tensor of the Vy-propagator contributes at leading order, since the second one is additionally suppressed 
as one can see from the defining condition ka + Ky < 2k av of the here considered case II. The relations between the 
IREs are 

5 A = -ka, (59) 
Sy = -kv, (60) 

$AV = K A V - K A - K V . (61) 

With this the new term in the formula for the lower bound on an IRE becomes 

28 av — 5a — <V — 2kav — 2ka — 2kv + ka + ny = 2kav — ka — «v > 0, (62) 

where we used the defining assumption of case II for the last step. Thus all diagrams that feature mixed propagators 
are additionally suppressed. This allows one to determine the value immediately, since if we disregard 

in the IR all diagrams that contain an yllZ-propagator or vertices that necessarily contain an ^4T^-propagator like the 
AAV-vertex, the usual Landau gauge system of DSEs remains. Taking into account that here the FT^-propagator 
has exactly the same form as the propagator of the fcrmionic ghost, the numerical factors from color and Lorentz 
indices add up in such a way that effectively only contributions from the Faddeev-Popov ghost contribute, as expected. 
Remember that we could not integrate out the y-ghost as the other ghosts, because it coupled to the gluon. In case II, 
however, this coupling is IR suppressed and does not play a role at leading IR order. The auxiliary field contributions 
cancel each other and the value of k is calculated as 0.595353 as for the standard Faddeev-Popov action. Since at 
leading IR order the DSEs are the same as in normal Landau gauge, we can use other results from it as well, e. g. the 
formula for the IRE of a vertex with 2n ghost- and m gluon- legs |38l 140] : 



K2 n ,m = (n - m)n + (1 - n) - 2j . (63) 

The IRE of the mixed two-point function can be calculated from its DSE. Here the importance of diagrams in the 
IR is determined by the number of AlZ-propagators, where one has to keep in mind that vertices exist that must have 
at least one mixed propagator in their DSEs. The one-loop diagrams of the I^A-DSE, as depicted in fig. [6] have the 
following IREs: 

Sa + Sav + k A aa, Sa+Sv + kaav, <V + <W + k A vv- (64) 

The fourth diagram existing in the yA-DSE is subleading, since it has two mixed propagators, kavv is zero due to 
the scaling relation and the IREs of kaaa and kaav can be determined from their DSEs as 

KAAA = 3(5y, (65) 

kaav = <W + 2<V- (66) 
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FIG. 6: Diagrams of the V^-DSE used explicitly in the text. 




Plugging this into eq. (64 1 and using eq. (59) as well as the scaling relation ka + % K v = 0, we get kav f° r all three 
diagrams, which proves that they are at leading IR order. For the calculation of kav we have to deal with several 
unknowns: The coefficients of the dressing functions (cIav, dv, cIa) and the IREs Ky and kav- However, cLav always 
drops out of the equations and dv and dA combine to dAdy, which can be calculated from the other propagator DSEs 
as 0.0267784. That leaves only kav ■ For a practical calculation we further truncate the system, since we cannot take 
into account all leading diagrams. We only take the AV -VV Aoorp into account, since it only contains AVV-vevt\ce&. 
The IRE of the AVV-veitex is 0, i. e. the bare diagram contributes to leading order. There are in principal other 
diagrams with the same order, but similar to the Faddeev-Popov ghost-gluon vertex, they can be neglected. They 
would correspond to a higher-loop correction. For the same reason we can neglect other diagrams in the y^4-DSE, 
since their IR leading contributions come from one-loop diagrams with IR enhanced vertices and plugging the IR 
leading parts of the vertices into the two-point DSE yields two-loop diagrams. In other words we take only into 
account the order dAdy and neglect everything of order (dAdy) 2 and higher. However, one should keep in mind that 
such truncations are only needed to calculate a numeric value for the IREs and that the qualitative features and the 
scaling relations themselves are derived for the full system without truncations. 

The equation to determine kav is obtained by projecting transversely in Lorentz space and by f abc in color space: 

(d -i)(P 2 )— dAV =- 9 *N [ ^((p+ g ) 2 )- i+ y)- i+ ;-( P v-(^) 2 ) ; (67) 



(27r) d V2p 2 d A d 



The coefficient dAV drops out of the equation. The final equation that has to be solved numerically is 
0.0513093r(l - K^)r(0.595353 + k A v) 0.0370825r(-K^)r(0.595353 + k av ) 



1 = 



r(1.40465 - K AV )T{1 + K AV ) r(1.40465 - k av )T(2 + k av ) 

0.0513093r(-K AV )r(1.59535 + k av ) 0.0210772r(-l - K A y)r(0.595353 + k av ) 



T(0.404647 - k av )Y(2 + k A v) r(1.40465 - K Ay )r(3 + k AV ) 

0.0247217r(-l - k av )T[1. 59535 + k A v) 0.0513093r(-l - k A v)T{2. 59535 + k A v) 

r(0.404647 - k A v)T{3 + k A v) r(-0.595353 - k A v)T(3 + k AV ) 

0.0237307r(-2 - K Ay )r(0.595353 + k A v) 0.0210772r(-2 - k av )T(1. 59535 + kav) 

r(1.40465 - K4v)r(4 + kav) r(0.404647 - K A y)r(4 + k av ) 

0.0370825r(-2 - K Ay )r(2. 59535 + kav) 0.0513093r(-2 - K AV )r(3. 59535 + k av ) 

r(-0.595353 - K4v)r(4 + k A v) + r(-l. 59535 - k^)F(4 + k av ) ' 



(68) 



As one can see from fig. [7] there are several solutions. A constraint for kav is ^av > Ky + ka = -Ky + d/2 — 2. 
The smallest solutions fulfilling this constraint are 0.0668776 and 0.981386. Further solutions always just a little bit 
below integer numbers follow. At this point we cannot say which one is correct. In our truncation all of them are 
admissible. 

An interesting observation is the fact that the AV^-FV-loop (see fig. [6]) is automatically at leading order: 

<V + &AV + KAVV = — H V + KAV — KA — Ky = KAV- (69) 



Therefore the tree-level term has to be canceled in order to allow a solution with k < 2. For the VV- two-point 
function we can employ the horizon condition in the same way as for the Faddeev-Popov ghost to explain the required 
cancelation of the tree-level term in the DSE, since their DSEs are the same in the IR. 
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FIG. 7: The plot to determine solutions for kav- On the Y-axis we plot the left- and right-hand sides of eq. (681 



C. Case III: ha + kv = 2kav, no cancelations 

Writing detC = c A c v + 2Nc 2 AV = (p 2 ) 2KAV detD = (p 2 ) KA+Kv det D, we get 

T-.AA,ab _ rab J_ p |__ / 7 rv\ 

~° p 2 ""(p2)«A det U' ^ 

nV^abcd _ ^ ^ /" ^acrhrf _ rabe rede p ^AV \ /71 \ 

~P 2 c y (p 2 )^ d 3a- / / ^"detDj' { ' 

n AV,abc _ ■ rabc J_ p ^AV , . 

J P * ^(p 2 )"^ detD' ^ 

If no cancelations occur cfet D is just a momentum independent constant and it assumes a similar role as the coefficients 
of the power laws, i. e. it is an unknown constant that has to be determined. The relations between the IREs are 

<5.4 = -ha, (73) 

S v =-k v , (74) 

5av = -K.AV- (75) 



The additional term in the formulas for the IREs, eqs. (26 1 and (27), becomes 

2Sav — 8a — Sv — —2kav + ka + kv = 0, (76) 

where we used the defining assumption of case III. Thus again the new terms vanish, i. e. we can directly get the 
scaling relation from the vertices in the Lagrangian, if we assume that the tree-level two-point function in the V^F-DSE 
gets canceled in a similar way as that for the other IR enhanced fields. The only non-trivial solution arises from the 
Afjr]- and AVV -vertices: 

d 

k:=kv = k v , ka + 2k=--2. kavv = K,A m — 0. (77) 

The IRE of the mixed two-point function can be calculated from the defining assumption as 

k d , 
k av = -- + --1. (78) 

In all DSEs the diagrams with a bare Af\r\- or j4VV-vertex are leading and the formula for the IRE of a vertex is 

K2n,m = (n - m)n + (1 - n) (- - 2 J , (79) 
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where m is the number ghion legs and n the number of legs of ghosts or auxiliary fields. 

Thus the qualitative behavior is the same as case II, i. e. the propagators of the Faddeev-Popov ghost and the 
auxiliary fields are IR enhanced and the gluon propagator is IR suppressed, as is the mixed propagator. The difference 
between case II and III is that the mixed propagator is more IR suppressed with respect to k in case II, since there 
it holds that Sav > k /2, whereas in case III it is k/2. 

Having found a scaling relation does, however, not yet mean that we can be sure there really exists a corresponding 
scaling solution. A first indication of its existence is a solution for the IRE k. This is an intricate task as the system 
of equalities that has to be solved involves the coefficients of the power laws, the dij, in a non-linear way. In case II 
only the IRE n appeared non-linearly and consequently one can obtain all solutions within a given interval. In the 
present case we were able to find several possible solutions with < k < 1, but all of them are in contradiction to 
requirements such as the positivity of d c and d A . Thus we cannot guarantee the existence of a solution in case III. 

For case III it is necessary that the tree-level term in the AV-DSE gets canceled by some mechanism that may be 
related to the horizon condition. If this is not the case the AV-VV-loop can only be leading for k > 2 or - if the 
tree-level term leads - we get k = 2. 



D. Case IV: ka + kv = 2k av, cancelations 

In the case that the IREs fulfill the relation K A + Ky = 2n A y it is possible that the leading terms in the determinant 
detC = c A c v + 2Nc 2 AV = (p 2 ) 2KAV detD = (p 2 ) KA+Kv det D cancel, i. e. d A d v = -2Nd Av . As the IRE of the next 
leading term is undetermined yet we introduce det D = du{p 2 ) KD , with kjj > 0. This bound comes from having 
factored out the IREs of the leading dressing functions in the determinant: The remaining pieces are constant or IR 
suppressed. In case IV the constant pieces cancel each other and something IR suppressed remains: 

detC = c A c v + 2Nc\ v = (p 2 ) KA+KV (d A d v + 2Nd\ v + d D (p 2 ) KD + . . .) = (p 2 ) KA+Kv (d D (p 2 ) K ° + . . .). (80) 

The propagators are then almost the same as in case III, except for the VV"-propagator, where now the first term 
is suppressed: 



n AA,ab _ cab 1 p dy 1 

» v p 2 ^(p 2 )^ detD 1 



T^VV,abcd 11/' xaccbd n rabe rede r> ^ AV \ ,1 1 / rabe rede D ^rf Av , 



T\AV,abc _ ■ fabc 1 p ^AV V * , R o\ 
~ S p 2 ^(p 2 )— detD' 



The relations between the IREs of propagators and two-point functions are in this case 



Sa = -ha - k d , (84) 
Sy = —Kv — k d, (85) 
<W = -K-AV - k d . (86) 



Thus we have again 2S A y — S A — 5y = 0. 



The inequalities relating the different IREs of propagators, eq. (18 1, and vertices, eq. (U7|, remain valid, since there 



does not appear any IRE of two-point functions. However, the appearance of kd influences the remaining power 



counting as follows. From eq. (27 1, which gives a bound from the analysis of propagator equations, we get for the 
VV^-propagator : 

K V + Sy - - (6 A + 2Sy) >0, 

k a + 2n v + k d >0. (87) 

We chose the yF-propagator because there is only one possible bare vertex in its DSE, contrary to the AA-DSE. On 
the other hand, from counting the VV-AA-loop in the VV-DSE, see fig. [5j we get another bound: 

Ky <Sy + S A + K A yy , 
K A + 2ny + 2k D <K A yy < 0, 
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where we used kavv < from the DSE of the AVV-vertex. Without kd we could now infer the usual scaling relation. 
However, combining these two inequalities yields now some additional information on K£>: 

—kd < ka + 2kv < —2kd, 

k d > 2k d - (89) 

This inequality cannot be fulfilled, since kd is non-negative and we are back at ko — (case III). Thus case IV does 
not yield a solution. 

V. CONCLUSIONS 

The analysis of the IR regime of the Gribov-Zwanziger Lagrangian as performed in this article requires a method to 
deal with the mixing of the fields at the two-point level. Some general formulas necessary for this, like the lower bound 
on the IREs of vertex functions, were given in Sec. |III| However, the complete analysis is severely complicated by the 
mixing of the fields, which leads to a matrix relation between two-point functions and propagators. This necessitates 
the distinction between four different cases, which assume different combinations of dressing functions to be more 
important than others. Only two of these assumptions allowed a consistent solution: We could explicitly exclude 
that the dressing functions of the mixed two-point functions dominate in the IR and thereby lead to a completely 
new IR behavior of the gluon and Faddeev-Popov ghost. The second dismissable case allowed for cancelations in the 
determinant of two-point functions. This would introduce a new IRE K£>, whose only possible value turned out to be 
zero, i. e. no cancelations take place. 

We found one solution that features a suppression of the contributions from mixed propagators compared to the 
other diagrams, i. e. all diagrams with mixing propagators can be neglected in the IR. This leads to a system of DSE 
that is in the IR exactly the same as in standard Landau gauge. As a consequence the calculated value for the IRE 
parameter n := Ky = k c is also the same, n ~ (93 - Vl201)/98 fa 0.595353. Using this result in the DSE of the mixed 
two-point function its IRE can be calculated. However, we did not find a unique solution but several allowed values. 
The lowest values are 0.0669 and 0.9814. For all values the AU-propagator is IR suppressed. An important feature of 
this solution is that the DSEs reduce to the same system as for standard Landau gauge. There is much information 
available on that solution that can be used here, e. g. the use of a bare Faddeev-Popov ghost-gluon vertex is known 
to be sufficient in numerical calculations with no immediate need to improve on the kinematic dependence. 

We found a second possible solution, where all IREs are strictly related by a scaling relation, k j 4 + 2k c = ka + 2kv = 
Ky + 2kav = 0, but the inherent structure of the equations did not allow us to find a value for the IRE parameter k. 
Nevertheless we can infer from the non-negativity of k c that the qualitative behavior of this solution is the same as 
that of the first one. In this case the relations between the IREs of the two-point functions and the propagators are 
just the negative of each other, e. g. ka = —5a- Both solutions found correspond to the same qualitative IR behavior 
of the propagators: The gluon propagator is IR suppressed, the ghost and auxiliary field propagators are IR enhanced 
with the same IRE and the mixed propagator is IR suppressed. In case II this suppression is more pronounced. 

In obtaining our results we had to improve the usual power counting, which normally does not need to differ 
between the tensors of Green functions. Because of the mixing, however, the corresponding dressing functions can 
combine in different ways for different tensors. Thus it was necessary to specify a concrete basis. Choosing the tensors 
that appear in the Lagrangian leads already to a UU-propagator with two different tensors. The (possibly) different 
behavior of these two tensors was important in determining the solutions presented here. 

In the present article we only worked with the conventional Gribov-Zwanziger action [?]. As this requires that 
the horizon condition is fulfilled the propagators of the Faddeev-Popov ghost and the auxiliary fields uj and U are 
necessarily divergent [35] . In functional equations this corresponds to a possible choice for the boundary condition of 
the ghost two-point DSE [21 . A different choice leads to the decoupling solution and of course the question arises if 
this solution could be realized somehow within the Gribov-Zwanziger scenario. It is worth reminding one here that a 
refinement of the Gribov-Zwanziger theory |29| provides in fact a way to get a finite non-vanishing gluon propagator 
and a finite ghost dressing function by introducing condensates of the auxiliary fields. Thus, a direct connection 
between the refined Gribov-Zwanziger scenario and the corresponding boundary condition in functional equations is 
to be expected. Zwanziger's argument about cutting off the functional integral at the Gribov horizon is not invalidated 
by imposing the decoupling boundary condition. However, an interesting difference between the Gribov-Zwanziger 
and the refined Gribov-Zwanziger scenarios is that in the latter no longer the configurations directly at the horizon 
dominate but configurations within the Gribov region. This can be seen by calculating the ghost self-energy a(k 2 ) 
[2"9"] . The study of the refined Gribov-Zwanziger theory and of the related decoupling solution within the functional 
set up presented here is under investigation. 

To continue and complete this work it will be interesting to try a complete numerical solution of the equations. 
The information obtained here will definitely be important input for this, as it is now known which diagrams are IR 
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leading and which are the relevant tensors. Even without the full solution available we can conclude that taking into 
account the restriction to the Gribov region in the path integral does not affect the IR properties as obtained from 
functional equations: The scaling relation remains valid and the auxiliary fields behave as the Faddeev-Popov ghost. 
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Appendix A: Extension of power counting to mixed propagators 



We present in this appendix the calculations to obtain eqs. (26 I and (27). They are extensions of calculations given 



in ref. |41| to the case of mixed propagators. For the convenience of the reader we also give an example of how to get 
the inequalities (17) and (18 1. 



1. IRE for an arbitrary diagram and lower bound on the IRE 



Eq. (26 1 is the formula for the IRE of an arbitrary diagram. It can be obtained by counting all the IREs of 
propagators and dressed vertices: 

k v =l~ +^2ni(Si - 1) + ^2 n i 1 ...i r ( K h...i r +Ctl...ir)+ 

i vertices,r>3 

+ n h...i r c ii...i r -Cv (Al) 

vertices,r>3 

Hi are the number of internal propagators with IR exponents Si, whereas the numbers of vertices are rii 1 .,,i r . Super- 
scripts d and b stand for dressed and bare, respectively. In case none is given, we refer to both. The sums ~}2 vertices r > 3 
extend over all vertices with at most r legs, rrii is the number of external legs. 

Without mixed propagators it is possible to use topological relations to get rid of the internal propagators. Here, 
however, the dependence on the mixed propagators will remain. To avoid cumbersome notation we restrict ourselves 
now to the case that we only have A- and ^-fields. The inclusion of the fermionic ghost field goes along the lines as 
presented in ref. (41] . The first relation expresses the number of loops in terms of the numbers of propagators and 
vertices: 

1= ^ n i+ n AV + l- ^ n h-ir- ( A2 ) 

i—A.V vertices,r>3 

Furthermore we express the number of internal A- and ^-propagators by 

n * = \( X! k i 1 ' " Vn 'i-*r _ m i ~ n AV J , i = A,V, (A3) 
where fej 1- " v denotes the number of times the field fa appears in the vertex fa x . . . fa r . Plugging these expressions 



".n j + I - n {{ - ^ n h...i r J ^ + a. w (A w - 1) - 

vertices,r>3 



i vertices,r>3 
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+ n AV 
This leads to 



E 2 I E fci 1 "' Vn ii-V ~ m i- n AV ] - 1) + 
-A, V \ vertices, r> 3 J 

E "■! - i ',)■ E (> D ^ !E 

vertices, r>3 vertices, r> 3 i 

^- 2 ) l 1 -^ E m * -^E m »^+ E <-v (-f +«i!...H.+2- J 

\ i—A.V J i vertices, r>3 

E^( E ^-^.^(1+4-1))+ E »?. ^-"+2- r 

\ vertices, r> 3 / vertices, 1 ) 

((s AV -i)-l E to" 1 )! ■ ( A4 ) 



2 i jL^ ■ ■ • 2 ' "* ' / ^— ' ' •■•« V 2 ' " 2 

£= -A.V \ t;er£zces.r>3 / vertices ,r>3 



k„ = - - 2 



1 - i E m * ) ~^E mi5i + \ n Av{^AV -6a- <V)+ 

i / i 

E <..«, ( (i - 1) (»■ - 2 ) + + \ E C"' -0 + 

vertices, r>3 V i / 

E iJlr 1 )^^^^)- 

ices,r>3 \ i / 



vertices, r>3 



(A5) 



In ref. |4T] it was shown that the last two terms are non-negative. We present one explicit example below. This 
conclusion does not change for mixed propagators. So the lower bound on the IRE of a diagram v is 



Kv,max = ( t! ~ 2 ) ( 1_ ^^ m *) - ^ ^ ]mjdj + ^n AV (2S A v -8a- 8 v ). (A6) 



In contrast to the non-mixing case it depends on the internal propagators, but only on the mixed ones. The value 
of 25av — 8 A — 8y has to be considered for each case. It should be non-negative, because otherwise the IRE is 
unbounded from below, since there are diagrams with an arbitrary number of mixed propagators. They can be 
explicitly constructed by appropriately replacing dressed quantities by their respective DSEs. 



2. Non-negativity of combinations of IREs 

For the IR analysis done in this article we also relied on inequalities obtained from RGEs. These functional equations 
also form an infinitely large tower of equations and are derived from the effective average action, see for example the 
reviews [T2"l IT^] . The average effective T k depends on the momentum scale k, above which all fluctuations are 
integrated out. The limits k — ¥ and k — > oo correspond to the standard effective and classical actions, respectively. 
The dependence of the average effective action on the scale k is given by: 

kd k T k [<f>]=±kd k R$ J G iJ . (A7) 

where dk '■— d/dk and the indices i and j represent field type, momentum and internal indices like Lorentz or color 
indices. Summation/Integration over indices appearing twice is understood. The connected two-point function is 

Gii = -crr4 u (A8) 

and the quantity Rjj is a regulator that implements the momentum cutoff at the scale k. T*- [$] is the corresponding 
two-point function. 
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RGEs are derived by functional differentiation of eq. (A7l with respect to fields. In contrast to DSEs they only 
feature one-loop equations, all quantities are dressed and a regulator insertion appears. Forms of regulators especially 
suited for an IR analysis are discussed in ref. [39 . 



The tower of RGEs allows one to derive the specific group of inequalities given by 



St > 0. 



(A9) 



We only provide a simple example here and refer for a complete proof to ref. [41 J. The simplest case is that of a 
three-point function. Let us therefore examine for concreteness the AVV- vertex. Part of its RGE is depicted in fig. [8] 
Power counting yields 



k avv <3kavv + 25y + Sa, 
0<kavv + ^(2Sv + S a )- 



(AlO) 



A corresponding analysis of its DSE would yield a similar inequality but with different, less restrictive numerical 
coefficients. Similar equations can be obtained for all vertices, as can be shown by an iterative proof [41] , and the 
set of inequalities can be written down in closed form, see eq. (A9l. One should note that these inequalities do not 
depend on the IREs of two-point functions and thus do not suffer from complications due to the mixing. Furthermore, 
we do not get additional inequalities if the fields mix. 
Another set of inequalities, given by 



(AH) 



can be derived by taking into account that the IRE of a vertex that appears in the action is non-positive, as can be 
inferred from its DSE. In the present example this means 



< 



-kavv <^ ( 2( V 



5 A ) 



25 x 



S A > 0. 



(A12) 



One should note that the number of inequalities in the second group, eq. (A9l, corresponds to the number of 



interactions in the Lagrangian, while the first group, eq. (All), is infinitely large. Furthermore, the interplay between 
DSEs and RGEs is nicely exhibited in the derivations of eqs. (A9l and (All I. The former can only be derived from 



the RGEs, while the latter require information provided by the DSEs. 
analysis. 



Yet, both inequalities are required for the IR 




FIG. 8: One specific diagram in the RGE of the AVV- 
vertex. Internal lines represent dressed propagators, black 
blobs dressed vertices. The gray blob is a regulator. 
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3. Analysis of leading diagrams in two-point DSEs 



For an arbitrary diagram in a propagator DSE we have from eq. (A5l 



k» = ( ^ - ^1 - i ^ - i ^2 m i S i + ^n AV (25 A v -S A - $v)+ 



dressed all 
vertices vertices 



if it is assumed to be part of the leading terms. Here «i can be any two-point IRE and Sj is restricted to IREs of 
non-mixing propagators. Using m, = 2 in the case of propagators, this can be written as 

«i + ^ E™^' ~ \ n Av(2S A V -S A -S V )- X "ii-ir I (4 - 1 j ( r - 2 ) + •', X'''' ' (i ' I = X n i-i r K ii-ir- 

j all \ 3 / dressed 

vertices vertices 

(A14) 

We can use the lower bound for the IREs of the vertices on the right-hand side given by the maximally IR divergent 



solution, eq. (26 1, to get a new inequality: 



J all V 3 J 

vertices 

dressed \ V 3 J 3 

vertices ^ ^ 



(A15) 



Here indicates the number of times a mixed propagator is contained in the diagram that determines the IRE 

of the vertex 4> il . . .<pi r - Note that k A 1 y lr only is different from zero for vertices that necessarily contain an AV- 
propagator like the AAV-vertex. The right-hand side depends on dressed vertices only, indicated by the d superscript 
of n. On the other hand, the left-hand side sums over dressed and bare vertices, so that in total only the bare vertex 
remains in the sums over vertices: 



m + l^mj <.., r ((± lj (r 2) + \ ]T kf-'-Si j - \(25 AV -8 A - 5 V ) 



( _ \ 

n AV + 



dressed 
vertices 



> o. 



(A16) 

J2i k l i 1 '" lr = r was used here. 
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